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Relaxation processes in collisionless dynamics lead to peculiar behavior in systems with long- 
range interactions such as self-gravitating systems, non-neutral plasmas and wave-particle systems. 
These systems, adequately described by the Vlasov equation, present quasi-stationary states (QSS), 
i.e. long lasting intermediate stages of the dynamics that occur after a short significant evolution 
called "violent relaxation". The nature of the relaxation, in the absence of collisions, is not yet 
fully understood. We demonstrate in this article the occurrence of stretching and folding behavior 
in numerical simulations of the Vlasov equation, providing a plausible relaxation mechanism that 
brings the system from its initial condition into the QSS regime. Area-preserving discrete-time 
maps with a mean-field coupling term are found to display a similar behaviour in phase space as 
the Vlasov system. 



I. INTRODUCTION 

The evolution of collisionless systems poses an interest- 
ing challenge in kinetic theory. Indeed, from the lack of 
a collision term in the kinetic equation ruling the evolu- 
tion of the system arises the need for another relaxation 
mechanism. Unusual relaxation properties are found in 
many physical systems among which self-gravitating sys- 
tems, non-neutral plasmas and wave-particle interactions 
for instance. These systems fall into the categories of 
long-range interacting systems, a domain of physics that 
is the object of renewed interest [THS]. 

One phenomenon in particular has been evidenced nu- 
merically Pl and theoretically [51 [B] : In systems with 
long-range interactions, increasing the number of parti- 
cles N causes the system to evolve, with a time scale of 
order 7V°, towards an intermediate state that is not the 
one predicted by statistical mechanics and whose life- 
time increases as TV'', with 6 > 1. These states are 
called quasi-stationary states (QSS) and are equilibria 
of the continuum limit given by the Vlasov equation [4]. 
QSS have been observed in the Hamiltonian Mean-Field 
(HMF) model [7], the free-electron laser [8] and also in 
self-gravitating systems [9]. Experimental perspectives 
regarding these QSS have recently been proposed [TO^ . 

The fast evolution on a timescale is termed vio- 
lent relaxation, following the terminology of the studies 
on ID self-gravitating systems. Early works report that 
the evolution towards a stationary regime is not driven 
by two-body encounters and question the nature of the 
initial evolution of gravitating systems [11] [T^] . In light 
of these observations, Lynden-Bell presented in 1967 his 
statistical mechanics of violent relaxation in stellar sys- 
tems (13j . The aim of his theory is to explain the outcome 
of collisionless dynamics and the unusual energy distribu- 
tion found in galactic dynamics. Lynden-Bell's (LB) the- 
ory is a statistical theory that takes into account, among 
other properties of the Vlasov equation, its incompress- 
ible character in the single-particle phase space; accord- 



ingly, an ergodic-type hypothesis on the dynamics allows 
one to compute, via an entropy maximization, station- 
ary states of the Vlasov equation. LB's theory predicts 
the magnetization in the HMF model and provides the 
most likely QSS solution for this model, as of now jl4j . 
The work of Lynden-Bcll has been followed by numerous 
numerical studies on violent relaxation, see Refs. P"51 - [T5] 
for instance. LB's theory or other attempts to compute 
the outcome of violent relaxation can be found in more 
recent work on the self-gravitating sheet model [3], non- 
neutral plasmas [I^ and 2D self-gravitating systems [20] ■ 
Recently, new macroscopic observables have been pro- 
posed to measure the approach to equilibrium in the self- 
gravitating sheet model [51]. An extensive study by the 
same authors indicate the physical situations in which 
LB's theory provides a good prediction and the reasons 
it does not work in other situations [22) . 

Let us also mention the use of the Vlasov equation in 
the field of plasma physics [23^ that is probably the widest 
area of research making use of it. The understanding of 
collisionless relaxation is also of great interest in this field 
and is mostly known for the famous problem of Landau 
damping (see Ref. [24] for instance). 

The purpose of this article is to demonstrate via direct 
numerical simulation of the Vlasov equation that stretch- 
ing and folding structures occur in some regions of phase 
space. A measure of the consequent deformation of the 
fiuid allows us to characterize the short-time evolution in 
the Vlasov equation. We choose to take as the main ve- 
hicle of our study, the Hamiltonian Mean-Field (HMF) 
model [25], a system of globally coupled rotators mov- 
ing on a circle. The HMF model has been the object of 
many studies as a paradigmatic representative of long- 
range interacting systems and the properties of its asso- 
ciated Vlasov dynamics are well-known. In addition, in 
the spirit of demonstrating fundamental dynamical prop- 
erties, the similarity of the HMF model with the forced 
pendulum makes it a worthwhile example, as evidenced 
in Ref. 
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Besides, we also consider discrete-time versions of 
Vlasov dynamics, in which the phase-space distribution 
function evolves under the effect of an area-preserving 
map. Two such mean-field maps are investigated, which 
are based on Arnold's cat map [571 [21] and the standard 
map j29l|30]. 

The plan of the paper is the following. In Sectionjllj we 
introduce the Hamiltonian Mean-Field model, its associ- 
ated Vlasov equation and the phenomenology of QSS. 



Section III presents the numerical algorithm and the 
computation of the perimeter that is used to substantiate 
quantitatively our claims. Section |IV] gives the results of 
Vlasov simulations, in which stretching and folding occur 
in phase space. In Section |Vj mean- field maps are stud- 
ied in order to determine the conditions of stretching of 
fluid boundaries in phase space. 



II. COLLISIONLESS DYNAMICS AND THE 
HAMILTONIAN MEAN-FIELD MODEL 

The Hamiltonian Mean-Field (HMF) model was intro- 
duced as a simplified model to study collective effects and 
has become a paradigmatic model to study long-range in- 
teractions O [25] . It is composed of particles on a circle 
interacting via a cosine potential and described by the 
following Hamiltonian : 
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where 9i is the position in the interval [— tt; tt[ of the ith 
particle, Pi its momentum and N is the number of parti- 
cles. At equilibrium, the HMF model is characterized by 
a second-order phase transition, identified by the magne- 
tization : 
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TO = |m| is equal to zero above the critical energy Uc = 
= |, while TO > for energies below Uc- 
Starting from an out-of-equilibrium initial condition, 
the HMF model displays interesting dynamical phenom- 
ena : in addition to the energy, the initial value of to 
influences the regimes attained by the system [T3], a de- 
pendence which is not found at equilibrium. The value 
that m reaches is not the one corresponding to equilib- 
rium statistical mechanics. This phenomenon lasts for a 
time lapse that depends on the number N of particles 
considered as iV^-^ after which the system eventually re- 
laxes to equilibrium. This intermediate regime is called 
a QSS [71IH]. 

Let us now introduce the Vlasov equation describing 
the evolution of the distribution function f{0,p) in the 
HMF model and valid in the thermodynamic limit |231 
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where e is the periodic spatial coordinate, p is the mo- 
mentum, V is the potential, depending self-consistently 
on /. The time evolution of Eqs. Q conserves the energy 
U : 

u[f]{t)= [ de dp f{e,p;t) 



X (^y + ^ (1 - m4f] cose- my[f] sin0)j , (4) 

the normalization J dedp f{e,p;t) = 1 and the total mo- 
mentum J de dp f{e,p;t) p. 

In order to solve numerically Eqs. ([s]), we use the semi- 
Lagrangian method (see Ref. [35] , or Ref. [33] for an ap- 
plication to the HMF model) with cubic spline interpo- 
lation. Although numerical limitations come into play 
j34[ . Vlasov simulations have proven useful to study the 
HMF model in the thermodynamic limit |31j . The semi- 
Lagrangian method displays a relatively small amount of 
numerical dissipation and thus fits adequately the pur- 
pose of computing the perimeter of the fluid. 

Lynden-Bell's (LB) theory has proven so far successful 
for the HMF model [14], the free electron laser [8] and 
partly for the self-gravitating sheet model [HI [22] ■ The ef- 
fectiveness of Lynden-Bell's theory to describe the HMF 
model has been an important confirmation of the ade- 
quateness of the Vlasov equation to understand QSSs. 

It has been recently shown that Lynden-Bell's theory 
captures fundamental properties of the Vlasov dynam- 
ics and that phase space stirring already provides a fair 
amount of evolution. The authors of Ref. [3S] obtain good 
comparison between Lynden-Bell's theory and a model 
with no interaction at all between particles. The same 
authors propose an exact solution to the mean-field equi- 
libria in the Vlasov equation, fully demonstrating how 
stirring allows the dynamics to reach steady states in 
Vlasov dynamics |36j . Stirring is however insufficient to 
obtain the better agreement found in the HMF model 
when using LB's theory. To complete the dynamical pic- 
ture, another relaxation mechanism is needed. 



III. EVOLUTION OF THE PERIMETER OF 
THE FLUID 

In the context of 2D fluid dynamics, experimental re- 
sults can be obtained by direct visual inspection. In 
the kinetic context, the distribution function (DF) lies 
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in phase space {9,p). The technique that we propose is 
reminiscent of methods used in fluid dynamics. To ana- 
lyze the complex behavior of the DF, we track its bound- 
ary in the same way as tracking a dyed region in a fluid. 
Instead of a dye, phase space is filled with a step profile 
of value /q. Inside the step, / = /o while outside of the 
step / = 0. An example of such step profiles is known 
as the waterbag (WB) initial condition. The WB initial 
condition possesses the following advantages : a simple 
formulation of LB's theory, a nice interpretation in terms 
of a dye in phase space and a very widespread use in the 
literature. 

The perimeter Pf of the DF / is the length of the in- 
terface between the f = fo and / = regions of phase 
space, it is similar to the intermaterial area per unit vol- 
ume defined in Ref. [37^ . Pf provides a direct measure of 
the deformation of the DF. 

The WB is defined by a width of 2A9 and 2Ap in the 
two dimensions of phase space. The initial perimeter Pf 
of this "fiuid" is 4{A9 + Ap). The preservation of phase- 
space volume in Vlasov dynamics implies that the area 
of the fiuid remains constant. 

In the simple case of free-streaming, it is easy to com- 
pute the time evolution of Pf : 



Pf{t) = AAe + AAp\/l + i2 



(5) 



which is seen to be asymptotically linear in time. This 
linear behavior is expected to hold in the case of regular 
dynamics, as in the case of integrable systems. 

In more complex situations, we expect a deformation of 
the contour of the fluid, similar to what happens in fluid 
dynamics [37] . Indeed, 2D flows submitted to parametric 
forcing exhibit chaotic behavior, and analogies between 
Vlasov dynamics and 2D Euler equations are common 
|38j . If neighboring fluid elements experience an expo- 
nential separation of their positions, we expect that the 
perimeter is dominated by this contribution so that 
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where A is a positive real number. After the stretching 
mechanism has taken place, the fluid folds onto itself to 
accommodate for the increasing perimeter at flxed area. 
In fluid dynamics, the obtention of an optimal mixing, 
via stretching and folding of the fluid, is a ongoing topic 
of interest [571155]. 

The perimeter is computed numerically with the help 
of the routine CONREC [40]. This routine, initially used 
to display contour lines on a computer display, has been 
modified in order to compute the length of this contour 
and has been implemented within the Vlasov simulation 
program used in Ref. J33] , a modification that is necessary 
in order to avoid storing all time steps of the DF for 
post-processing. Equation ^ is checked to give a perfect 
match with the simulation for m — (data not shown). 
A check of our conjecture for non-interacting systems is 
tested for a pendulum (a system similar to Eqs. (|3| but in 
which the magnetization m is set to a constant value in 



the equations of motion). Figure [I] displays the evolution 
of the DF and Fig. [2] displays the evolution of Pf (t) for 
the same simulation. 




FIG. 1: Vlasov simulation for a non-interacting set of pen- 
dula. The top panel displays the waterbag initial condition, 
parametrized by — 2.85 and Ap = 0.79. The bottom 
panel displays only the contour at time t = 10. The filamen- 
tary structure is caused, in this system, only by the differential 
rotation in phase space, i.e., phase-space stirring. 
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FIG. 2: Perimeter of the fluid as a function of time for the 
simulation shown in Fig. [T] The top panel uses a logarithmic 
scale on the j/-axis and the bottom panel uses a linear scale 
for the 3/-axis. 



Table |T] lists the numerical parameters for the Vlasov 
simulations. All initial conditions are runs with all these 
parameters, allowing to monitor the convergence of the 
physical quantities as a function of the number of grid 
points. 
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TABLE I: Numerical parameters for the Vlasov simulations. 





Ng 




At 


0k5 


512 


512 


0.1 


Ik 


1024 


1024 


0.1 


2k 


2048 


2048 


0.1 


4k 


4096 4096 


0.1 


8k 


8192 


8192 


0.1 


8kb 


8192 


16384 0.1 



IV. STRETCHING AND FOLDING IN PHASE 
SPACE 

The present section focuses on the interacting HMF 
model. The variations of the self-consistent field m cause 
a complex evolution of the initial waterbag that we con- 
sider in several test-cases. The first case (run f in Ta- 
ble |ll| leads to a magnetized state and has the property 
that keeps its sign in the course of time. Also, a 
significant amount of mass remains enclosed inside the 
separatrix at all times, meaning that it is not subject to 
alternating stretching motion but only to a stirring mo- 
tion that is reminiscent of the non-interacting situation. 
A variation on run f is the run j that also leads to a mag- 
netized regime and possesses a similar structure in phase 
space, mo for run j is however equal to the predicted 
value of TO in the QSS regime by Lynden-Bell's theory 
and we expect a smaller amount of stretching and folding 
because the variations of to are smaller. The second case 
(run i in Table refers to the peculiar value U — 0.69, 
with Too = 0.20, which is known to possess counter ro- 
tating resonances, then experiences changes of sign 
that imply elliptic-hyperbolic bifurcations [38] which we 
discuss further on in section HVBl 

TABLE IL Parameters for the waterbag initial condition of 
runs f , i and j . 



Run 


U 


mo 


f 


0.60 


0.20 


i 


0.69 


0.20 


j 


0.565 


0.50 



The organization of phase space is similar to the one 
of a pendulum at every given time. The simplest situa- 
tion is when niy = at all times with small variations 
of TTix around a finite average value. The system is then 
close to a set of pendula, although small oscillations of 
rrix are sufficient to induce stretching and folding. In 
the presence of stronger oscillations of rrix , the observed 
phenomenon of stretching and folding is increased. The 
oscillations of the separatrix, which drive the stretching 
and folding, are indeed linked to the one of m^. The re- 
gion of phase space that is located under the separatrix at 



all times experiences deformations, but the organization 
of phase space is such that around the elliptic fixed point 
at (0, 0) a part of the waterbag remains intact, reminding 
the "core-halo" phenomenon observed in self-gravitating 
systems PE2]. 

In the counter propagating resonances regime, first ob- 
served in Ref. [25' , an example of hyperbolic-elliptic bi- 
furcation is found. Such a behavior has been observed 
already in the single- wave model in Ref. (SSj. The au- 
thors of the latter reference find that a succession of 
elliptic-hyperbolic bifurcation produces strongly chaotic 
Lagrangian trajectories and destroys dipolar structures, 
modifying the general behavior of the system. As 
it stands for the HMF model, two coherent counter- 
propagating clusters remain present in the system. 

A. Variation of tUx at null rriy 

We consider run f of Table [ll] in this situation, my — 
at all times and rrix displays strong initial oscillations 
from t — to t Ri 100 followed by smaller oscilla- 
tions (see Fig. [s]) around a finite value. The initial time 
lapse corresponds to violent relaxation, during which we 
are interested in the behavior of the fluid. The initial 
stronger oscillations are typical of violent relaxation in 
self-gravitating systems with the difference that the mon- 
itored quantity in these systems is the virial ratio [15j . 
The value of m, in the HMF model, is sufficient to fully 
follow macroscopic quantities. 

0.50| < < < 1 




FIG. 3: Evolution of m{t) as a function of the time t for run f 
of Table [ll] for different numerical settings (see Table |l]l . The 
initial behavior is very similar for all runs, later times display 
small phase and amplitude differences. 

Figure [s] displays m(t) as a function of time, for run 
f and the parameters of Table |Tj Between t = and 
t « 50, the curves superimpose well. Up to t fa 100, 
qualitative agreement is still valid. We may already con- 
clude that the amplitude of the oscillations decay, a key 
point in the theory of violent relaxation. However, the 
amplitude of the oscillations in m(t) behave in a different 
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manner than found in Ref. [33] where increasing the num- 
ber of grid points enhance sustained oscihations. Here, 
increasing the number of grid points does not lead to a 
monotonous increase in the amphtude of the oscillations. 
We can quantify this behavior by taking the standard 
deviation of the time series m{t), taken here between 
t = 100 and i = 200. The resuhs are given in Table [111] 
The standard deviation is found to decrease from param- 
eters "0k5" to "4k" and to increase afterward. This can 
be understood by the presence of two competing reasons 
for the damping. A low number of grid points increases 
the phenomenon of numerical dissipation (as found in 
Ref. j33j) while a high number of grid points enables the 
description of the stretching and folding in phase space. 
This latter phenomena, by allowing a wider spreading of 
fluid elements, contributes to the homogeneity of the sys- 
tem, from the point a view of a given energy level, and 
thus to a reduction of the oscillations in m{t). 



TABLE III: Standard deviations of m{t) for run f with differ- 
ent parameters. m{t) is taken between t = 100 and t = 200. 



0k5 


Ik 


2k 


4k 


8k 


8kb 


2.5 10"^ 


1.7 10"^ 


1.3 10" 


1.1 10"^ 


1.5 10~^ 


1.6 10"^ 



A general view of phase space is given in Fig. |4] The 
phase-space structure is as follows : there is an elliptic 
fixed point at {0,p) — (0, 0) and an hyperbolic fixed point 
at {9,p) = (±7r, 0). A separatrix joins the two hyperbolic 
points, but its height is variable as the magnetization 
evolves in time. This time dependence allows one to con- 
sider the HMF model as a time-dependent pendulum, 
keeping in mind that in the HMF model the particles 
generate self-consistently the value of m. Three differ- 
ent behaviors are found, depending on the energy with 
respect to the separatrix, consistent with the findings of 
the forced pendulum [35] : 

• Below the separatrix energy, the so-called trapped 
particles remain close to the elliptical point and 
perform an oscillatory motion. 

• Above the separatrix, particles cross space with a 
velocity of constant sign. 

• At energies near the separatrix, particles experi- 
ence a very complex behavior, known as separatrix 
crossing. 

The oscillations of the separatrix will cause fluid elements 
to experience stretching and folding. 

Figure [s] displays Pj (t) for the parameters from Ta- 
ble |T] Increasing the number of grid points allow to reach 
a higher value of Pf before saturation. Before satura- 
tion, all runs show quantitative agreement on the value 
of Pf{t), indicating convergence as the number of grid 
points is increased. The initial behavior of Pf(t) is ex- 
ponential (notice the logarithmic scale of Pf in Fig. [5]), 
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FIG. 4: Phase space for run f at time t = 25, with the setting 
"8k", the contour is displayed, the region in red corresponds 
to / = /o and the region in white corresponds to / = 0. The 
instantaneous separatrix is drawn and three separate regions 
are visible : well inside the separatrix, around the separatrix 
and well above the separatrix. The filamentary structure is 
the result of the self-consistent evolution of the system. The 
box indicates the zoomed region of Fig. [6] 

showing that part of the perimeter is experiencing an 
early exponential separation. 




t 

FIG. 5: The perimeter Pf{t) for run f, for several configura- 
tions of Tablejl] The initial behavior is an exponential increase 
of Pf as a function of time (notice the logarithmic scale for the 
y-axis), followed by a saturation due to a numerical limita- 
tion. It is observed that increasing the number of grid points 
increases the saturation level. 

We notice that the stretching rate A of the perime- 
ter differs from a Lyapunov exponent in several respects. 
The perimeter is associated with a certain distribution 
function so that its stretching rate may depend on the 
choice of the initial distribution function. The perime- 
ter may undergo a regime of exponential growth, which 
does not go on, so that the stretching rate may not be 
defined in the long-time limit. Nevertheless, the concept 
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of Lyapunov exponent characterizing chaotic dynamics is 
useful to interpret the present stretching rate. In Fig. [5] 
the early exponential increases observed for Pf can be 
understood in terms of the phase-space dynamics of the 
pendulum periodically driven by the large early oscilla- 
tions of the mean field seen in Fig. [3] Indeed, the peri- 
odically driven pendulum is known to be chaotic, which 
can induce the stretching of phase-space domains leading 
to a positive Lyapunov exponent. Therefore, the perime- 
ter can undergo an exponential growth as long as this 
induced stretching goes on. However, as time increases, 
the oscillations of the mean field become of lower ampli- 
tudes as observed in Fig.jSj which reduces the extension 
of the chaotic zones in the phase space of the periodically 
driven pendulum and tends to decelerate the growth of 
Pf as seen in Fig.[5j 

Finally, to illustrate the behavior at small scales, we 
display in Fig. [6] a zoom of Fig. [1] at a later time. The 
stretching and folding structures are clearly apparent. 
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FIG. 7: Phase space DF f{9,p) for run i at time t = 20 
and t = 25 (see Fig.[8|. Two counter-propagating clusters are 
found that lead to sign changes in rrix . The top panel displays 
a separatrix for rux < and the bottom panel displays a 
separatrix for > 0. 
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FIG. 6: Zoom on Fig.|4]at a later time t = 100. The complex 
evolution of the initial waterbag leads to a stretching and 
folding behavior. A very intricate sequence of / = and 
f = fo regions can be observed. 




FIG. 8: nix and m as a function of time for run i. The vertical 
lines indicate the times at which the snapshots of Fig. [7] are 
taken. 



B. Elliptic-hyperbolic bifurcation 

The study of run i, leading to the counter propagat- 
ing resonances, displays an interesting bifurcation of the 
elliptic fixed point (0,0). The nature of that point de- 
pends on the sign of rrix that is found to change during 
time evolution. The separatrices corresponding to the 
two situations of the elliptic-hyperbolic bifurcation are 
shown in Fig. [7j 

The evolution of m and nix are displayed in Fig. [8] 
{niy = at all times). Fluid elements that are located 
in the alternating separatrices region experience a strong 
amount of stretching and folding. This is confirmed in 
Fig. [9] where Pf (t) is shown and compared to the result 
for run f. The improved mixing in run i should not how- 
ever lead to premature conclusions about the relaxation 



of the systems on a macroscopic level. The counter prop- 
agating resonances are very stable and prevent the re- 
laxation to complete equilibrium in which the resonances 
are not predicted. 

In order to link the phenomenon of stretching and fold- 
ing to the amount of evolution from an initial condition 
to a QSS, run j provides an interesting comparison. The 
initial magnetization is mg = 0.5 and the magnetization 
predicted by Lynden-Bell's theory is also mgss — 0.5. 
A large region of the initial waterbag will remain below 
the separatrix and smaller oscillations than for run f are 
expected. This is indeed found to be the case in Fig. |9] 
A comparison between run f and run j for the numeri- 
cal parameters 8kb is given in Fig. [9] confirming that the 
amount of stretching and folding needed for run j is well 
below the one for run f. 
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t 

FIG. 9: The perimeter Pf{t) for the runs f, i and j of Table |Tl] 
All runs are performed with the setting "8kb" of Table [I] 



A. Mean-field cat map 

In this case, the map takes the following form: 
Pn+i ^ Pn + Kxn - TOj; [/] sin(27ra;„) + my[/] cos(27ra;„) 

Xn+l Xn + Pn+1 (modulo 1) , 

(9) 

with the mean field 
m[/] = [/] +tmy[f]^C J dx dp f{x, p) e'^-- . (10) 

Moreover, we suppose that the initial distribution func- 
tion f{x,p) is uniform in a circle of radius j centered on 
the point x ^ p = | of the unit square, which defines the 
initial domain Dq. After n iterations, the distribution 
function is still uniform but in a domain D„ resulting 
from the time evolution of the map Q. 



V. MEAN-FIELD MAPS 

In this section, we simplify the dynamics of the HMF 
model into discrete-time maps giving the mean-field ap- 
proximation of symplectic coupled map systems |41j 
ruled by the periodically kicked Hamiltonian: 
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In the mean-field approximation, the positions and mo- 
menta of the particles after each kick are mapped by the 
time evolution according to 



Pn+l =Pn - U^iXn) - J dxdpf{x,p) U' {x^ - x) 
Xn+l = Xn + Pn+1 , 



(8) 

where n denotes the discrete time and f{x,p) the dis- 
tribution function representing the ensemble of particles 
normalized by J f{x,p)dxdp — 1. The map is area- 
preserving so that the normalization condition is main- 
tained during the time evolution. Here, we consider the 
interaction U{x) — C [1 — cos(27ra:)] and the external po- 
tentials: Uq — Kx"^ /2 corresponding to Arnold's cat map 
dZl US] or C/o = A [1 - cos(27ra;)] /(27r) giving the stan- 
dard map EaEO]. 

We suppose that the distribution function f{x,p) is 
uniform in a phase-space domain £)„, which evolves in 
time under the mean- field mapping D„ — ^"^{Dq). The 
area of this domain is preserved: Area(D„) = Area(Do)- 
The problem is now to determine the time evolution of 
the perimeter of the domain £)„. 



FIG. 10: The perimeter Pf{t) for the mean-field cat map 
with K = 0,1,2 and C = 10. The calculations are performed 
with 10^ (empty symbols) and 4 x 10^ (full symbols) points 
to represent the boundary and the bulk of the domain where 
the distribution function is uniform. 

If the coupling parameter vanishes, C — 0, the mean 
field disappears and the map reduces to Arnold's cat 
map, which is known to be fully chaotic if the param- 
eter X is a positive integer [27l [28] . We notice that the 
hyperbolic character of the map persists as long as the 
coupling parameter C is small enough that K ^ 27r|m|. 
Under this condition, we may expect that the hyperbolic 
character of the map will tend to stretch the domain D„ 
into a filamentary phase-space structure, which produces 
an effective uniform distribution over the unit interval x 
such that the mean field (10 1 tends to vanish as n — > oo. 
For the same reason, the perimeter of the domain Z?„ 
grows exponentially in time at a rate close to the Lya- 
punov exponent of the Arnold cat map. This effect is 
indeed observed in Fig.[lO| which depicts the perimeter 
of the domain _D„ versus the number n of iterations for 
the parameter values C = 10 and K = 0,1, 2. The com- 
putation is performed for two approximations in which 
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the boundary of the domain D„ and its interior are dis- 
cretized into 10^ and 4 x 10^ points. In the hyperbohc 
regime for i^T = 1 and K = 2, the growth of the perime- 
ter is exponential. However, in the absence of chaos for 
K = 0, the perimeter no longer increases exponentially 
and the computation is more sensitive to the discretiza- 
tion. The fact is that for the value C = 10 of the coupling 
parameter, the mean field vanishes so that an exponen- 
tial stretching of the perimeter is not even induced by 
the mean field. 



B. Mean-field standard map 

Now, we consider the mean-field standard map defined 
on the unit square as 



non-vanishing for < C < Ct by a phenomenon ana- 
logue to the one happening in the HMF model. For 
C = 1, the mean field fluctuates around the mean value 
nix — —0.21 with my = 0. This explains the early ex- 



ponential growth of the perimeter observed in Fig. 11 for 
A = and C = 1. Indeed, if we supposed that the mean 
field remains essentially constant, the Vlasov dynamics 
would be equivalent to the one of the simple standard 
map for A = —nix ~ 0.21. In this case, the early growth 
rate would be about A ~ 0.44, which is close to the value 
A ~ 0.48 of the mean-field dynamics with A = and 
C = 1 shown in Fig.[TT| 

These results confirm that the stretching of the perime- 
ter by the Vlasov dynamics is controlled by the local 
hyperbolicity of the effective single-particle phase-space 
dynamics induced by the current value of the mean field. 



■ Pn 



rriylf] cos(27ra;„) 

^ (11) 
with the same mean field as in Eq. ( 10 ). Here, the single- 



Pn+l 
Xn+1 



(A 



iTix [/])sin(27ra;„) 
(modulo 1) , 



particle potential as well as the one induced by the mean 
field have a similar trigonometric form so that the mean 
field TUx appears to have an effect comparable to the one 
of the constant parameter A of the standard map (in 
the case where my = 0) [3^1 [3U]. The initial distribution 
function f(x,p) is the same as in the previous subsection. 

Figure [TT] shows the growth of the perimeter of the 
domain _D„ as a function of the discrete time n for the 
parameter values C = 1 and ^4 = 0, 1, 2 for a discretiza- 
tion of the boundary and bulk of the domain Z3„ into 10^ 
and 4 x 10^ points. The early growth of the perimeter un- 
til about n = 6 iterations is well approximated by these 
discretizations for A = 1,2 and until about n = 15 iter- 
ations for A = 0. Thereafter, both approximations devi- 
ate, showing the sensitivity of the calculation to the dis- 
cretization of the domain Z3„ and the distribution func- 
tion f{x,p) into points. Nevertheless, the early growth 
appears to be exponential for A — 0,1, 2. 

For A = 1,2, the exponential increase of the perimeter 
can be explained by the global chaoticity of the stan- 
dard map for A > Ac 0.15 [HI ED]- Indeed, the stan- 
dard map is stretching domains at a rate equal to the 
Lyapunov exponent A ~ In^nA), which is positive for 
A = 1,2. At each iteration, the perimeter of the domain 
Dn is stretched by a factor exp(A) > 1 for every point 
in the chaotic zones of the map. The domain I?„ thus 
develops a filamentary phase-space structure uniformly 
distributed over the unit square so that the mean field 



( 10 1 tends to vanish and the mean-field standard map 



(111 behaves on long time as the single-particle standard 
map with A — 1,2 and m = 0. Accordingly, the expo- 
nential growth of the perimeter of the domain £)„ could 
go on as the discrete time n increases. 

In contrast, for A — 0, the dynamics is entirely de- 
termined by the mean field (10 1. It turns out that the 
mean field vanishes for C > Ct with Ct ~ 2, but is 




D 5 10 15 20 

iteration 

FIG. 11: The perimeter Pf{t) for the mean-field standard map 
with K = 0,1,2 and C = 1. The calculations are performed 
with 10^ (empty symbols) and 4 x 10*' (full symbols) points 
to represent the boundary and the bulk of the domain where 
the distribution function is uniform. 



VI. CONCLUSION 

Measuring the degree of mixing in kinetic theory is 
undoubtedly a considerable challenge. However, in order 
to assess the evolution of a given system, a quantitative 
tool must be chosen. An analogy with fluid mechanics 
allowed us to consider the use of the perimeter of the 
fluid in phase space, viewed as a dyed region. 

In both the Haniiltonian Mean-Field model and the 
mean-field cat and standard maps, our numerical exper- 
iments report an exponential growth of the perimeter, 
sign that exponential stretching occurs at least locally 
on the contour of the fluid. The fluid in phase space dis- 
plays folding and stretching structures that are absent 
from similar non-interacting and non-chaotic systems. 

The complex intertwining of / = and / = /o regions 
in phase space conflrms that a coarse grained theory is 
able to provide a good description of the fluid after mix- 
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ing has taken place. The existence of an infinite number 
of conserved quantities in the Vlasov equation forbids 
the increase of entropic functionals, which imphes that 
coarse graining is needed in order to describe the evo- 
lution towards stationary states (SS) or quasi-stationary 
states (QSS). In the light of our findings about the mech- 
anism leading to such QSS, violent relaxation can be un- 
derstood from a pure Vlasovian picture in which the de- 
formation of the initial condition allows the system to 
reach a QSS. This deformation takes place at constant 
phase-space volume but with a perimeter that increases 
exponentially during some lapses of the time evolution. 
The increase in the perimeter is made possible through 
successive stretching and folding the fluid. Subsequently, 
finite N effects may come into play but are of a different 
type and are not discussed herein. 

The understanding of phenomena related to violent re- 
laxation, for instance the separation of phase space in dif- 
ferent regions ^42] that leads to the formation of a core 
and a halo are coherent with the results of the present 
article. Applying our method to more common systems, 
such as the self-gravitating sheet model, is expected to 
reveal a similar behavior and would prove interesting to 
complete the existing body of literature. 

From our numerical observations, we have also found 
that stretching and folding cause a damping of oscilla- 
tions (see Table III) that competes with numerical oscil- 
lations. 

As known from the literature, the presence of self- 
consistent resonances may organize phase space in a 



way that prevents the statistical prediction of stationary 
regimes. Indeed, these resonances imply that the energy 
distribution function is non monotonous, in contrast to 
the predictions of Lynden-Bell's theory or Boltzmann- 
Gibbs equilibrium. The self-consistent resonances are of- 
ten at the origin of the so-called "dynamical effects" that 
prevent a system to reach a QSS or thermodynamical 
equilibrium. 

Finally, let us recall that, in the perspective given in 
the present paper, the phase-space exponential separa- 
tion leading to mixing has similarities to what happens in 
chaotic dynamics although notable differences exist such 
that the dependence on the initial distribution function 
as well as the time dependence of the stretching rate. A 
possible way to tighten both aspects is to use test par- 
ticles that feel the time-dependent magnetization gener- 
ated by a Vlasov simulation. It is possible to push this 
method further by the use of a Poincare map on these 
test particles, providing one with the tools of iterative 
mapping, which are of common use in nonlinear dynam- 
ics. 
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